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Abstract. Mutually unbiased bases (MUBs) are considered within the framework 
of a generic star-product scheme. We rederive that a full set of MUBs is adequate 
for a spin tomography, i.e. knowledge of all probabilities to find a system in 
each MUB-state is enough for a state reconstruction. Extending the ideas of 
the tomographic-probability representation and the star-product scheme to MUB- 
tomography, dequantizcr and quantizer operators for MUB-symbols of spin states and 
operators are introduced, ordinary and dual star-product kernels are found. Since 
MUB-projectors are to obey specific rules of the star-product scheme, we reveal the 
Lie algebraic structure of MUB-projectors and derive new relations on triple- and 
four-products of MUB-projectors. Example of qubits is considered in detail. MUB- 
tomography by means of Stern-Gerlach apparatus is discussed. 
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1. Introduction 

Since the early days of quantum mechanics, much attention has been paid to a 
problem of a good description of quantum states. The notions of wave function if) 
and density matrix p are most widely known and used. Nevertheless, these notions 
give rise to a problem of interpretation, especially, in case of measuring a quantum 
state. Outcomes of quantum observables are known to be probabilistic. In view of this, 
quasiprobability distribution functions like Wigner PU-function [1], Sudarshan-Glauber 
P-function [21 [3], and Husimi Q-function [3| are often used in quantum optics along 
with the wave function and density matrix formalism. The main drawback of W, P, 
and Q functions is that they cannot be measured experimentally. The problem of 
measuring quantum states resulted in developing quantum tomography, and then in a 
tomographic-probability representation of quantum mechanics (historical background is 
given in the review [5]). According to such a representation, any quantum state of light 
is described by measurable tomograms: optical, symplectic, and photon-number ones 
(see, e.g. the review [6]). As far as a finite dimensional Hilbert space is concerned, one 
can alternatively utilize spin tomogram [3, [8] and spin tomogram with a finite number 
of rotations [9]. Apart from being appropriate for reconstructing the density matrix, 
quantum tomograms are themselves notions of quantum states. Within the framework 
of the tomographic-probability representation, operators are described by tomographic 
symbols satisfying rules of the corresponding star-product scheme. We cannot help 
mentioning some of the probability-based approaches to quantum mechanics, namely, 
the expectation-value representation [TDl EH] and the Bayesian interpretation [T21 [13] 
utilizing a symmetric informationally complete positive operator-valued measure (SIC- 
POVMs are discussed, e.g., in pU UM HE]). 

The aim of this paper is to develop the star-product [17] quantization scheme based 
on mutually unbiased bases (MUBs) [T8l [19]. MUBs represent themselves a highly 
symmetrical structure and have many interesting properties (see, e.g. [20] and references 
therein). For example, a full set of MUBs is known to exist whenever the dimension 
of Hilbert space is a prime number or the power of a prime. We consider neither the 
problem of existence of MUBs in a given Hilbert space nor the problem how many MUBs 
there exist. We assume that the full set of MUBs is known for the space involved. In this 
article, we combine MUBs with the tomographic-probability representation. As a result, 
MUB-tomography of spin states is introduced, MUB-symbols of quantum operators are 
considered within the framework of the star-product scheme, the Lie algebraic structure 
of MUBs is pointed out, and new properties on MUB-projectors are derived. Special 
attention is focused on qubits. 

The paper is organized as follows. 

In Section [2} MUBs are briefly reviewed. In Section [31 MUB-based tomography 
is considered, scanning and reconstruction procedures are presented. In Section [U we 
follow the ideas of a generic star-product scheme [2H [22] and analyze the star product of 
MUB-symbols. In Section [5], an example of qubits is considered in detail. In Section [61 
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a practical realization of MUB-tomography by means of Stern-Gerlach apparatus is 
discussed. In Section conclusions and prospects are presented. 

2. Mutually unbiased bases 

Let us consider a d- dimensional Hilbert space endowed with a full set of mutually 
unbiased bases. If this is the case, MUBs consist of d + 1 bases {|act)}a=o ; where 
a = 0, . . . , d is responsible for the basis number, an index a = 0, . . . , d — 1 refers to one 
of the basis states belonging to the particular basis a. MUBs are to satisfy the following 
property: 

2 1 

\(aa\b/3)\ = -(1 - 8 a)b ) + 5 a>b 5 a ^, (1) 

where 5 a ,b is a Kronecker delta symbol. Eq. ([T]) implies that each basis is orthonormal 
and arbitrary two states belonging to different bases are equiangular, i.e. \{aa\b[5)\ 2 = k 
if a 7^ b. 

Let us now consider rank-1 MUB-projectors Il aa = \aa){aa\. 

Obviously, operators fl aa are semi-positive and satisfy the trace relation 
Tr[n aQ ,n b/3 ] = i(l — <5 a 6 ) + 5 aj bO~ a ,i3- An immediate consequence of orthonormality is 

d-l 

n a a = / for all a = 0, . . . , d, (2) 

d d-l 

^^n aQ = (d + i)/, (3) 

a=0 «=0 

where J is the identity operator. Since the relation ()2]) is valid for all a = 0, . . . , d, the 
total number of linearly independent operators Il aa equals 1 + (d+ l)(d — 1) = d 2 . This 
means that the identity operator / together with d 2 — 1 operators {LI aQ }, a — 0, ... ,d, 
a = 0, . . . , d — 2 form a basis in <i-dimensional Hilbert space. As a result, any operator 
including the density operator p of a quantum state can be resolved through these basis 
operators. Indeed, 

d d-2 

p = cjl + ^ XI ^/^/J, ( 4 ) 

6=0 /3=0 

where cj and are real parameters because both p and fl^ are Hermitian. Applying 
the trace operation to both sides of (BJ and utilizing Trp = TrLT^ = 1, we readily obtain 




«=-i i-v>:<*s . (5) 



>=3'+E5>(a*-3')- w 

6=0 /3=0 v 7 

Taking into account the non- negativity of operators Ii aa and the sum rule ([3]) , it is not 
hard to see that operators E aa = (c/+l) _1 II aQ , altogether form a positive operator- valued 
measure (POVM). We will refer to such a POVM as MUB-POVM. 
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3. MUB-tomography 



In an experiment, probabilities of measurement outcomes are only accessible. 
Tomography is a procedure allowing one to reconstruct density operator p with the help 
of measured probabilities. We will consider projective (von Neumann) measurements 
associated with MUBs. In other words, we assume that the probabilities 

p aa = (aa\p\aa) = Tr[pfl aQ ] (7) 

are known for all a = 0, . . . , d, a = 0, . . . , d — 1. As a consequence of expressions ©-Q 
we obtain the following normalization conditions: 

d-1 

^ ] Paa = ■ 



a=0 



d d-l 

EE 

a=0 a=0 



Pa 



d+l. 



(8) 



The physical meaning of p aa is the probability to find a system in the state \aa) 
which is itself an element of MUBs. The problem is to express the density operator p 
through probabilities p aa . This problem is solved in [23]. For the sake of the subsequent 
consideration we rederive the result and present it in a slightly different manner. 

Proposition. A reconstruction procedure of MUB-tomography reads 

d d-l 



P 



1 



d+l 



J2 J2 p w 

6=0 0=0 ^ 

Proof. Multiplying both sides of Eq 
we obtain 

1 

! d 



(9) 



§6§ by Ii aa and applying the trace operation, 



Pa 



where M n 



d d-1 

EE 

6=0 0=0 



(10) 



L a a,b/3 = Tr [11^11^] — \ = ^a,b(^a,i3 — ^) • The composed index aa can be 
considered as a single one: k = ad + a, k = 0,...,d 2 + d — 1. Then formula (fTUj) is 
nothing else as a linear system of equations with respect to c^, and can be rewritten in 
matrix form as follows: 



/ po - l/d \ 




( c 


\ 




= M 






\ PaV+d-1 ~ l/d ) 




\ CcP+d-l 


J 



where M is a {d 2 - 1) x (d 2 - I) 



M 



( M 






M 



block-diagonal matrix of the form 
\ 







v 







(12) 



\M ) 

with each {d — 1) x (d — 1) block A4 being equal to 

/ 1 - l/d -l/d ■ ■ ■ -l/d \ 
-l/d 1-1/d ■■■ -l/d 



M 



\ -l/d -l/d 



(13) 



1-1/d ) 
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A direct calculation of the inverse matrix yields 

/ 2 1 -•• 1 \ 
1 2 ■■• 1 



(14) 



\ 1 1 ... 2 J 

from which it is not hard to find an explicit solution of system ( flOj) . The result is 



Cb/3 = Pbf3 



n d-2 



13=0 



Pb(3 ~Pb,d-l- 



(15) 



The right-hand equality in (1151) is due to normalization condition (jHj). Substituting the 
obtained value of Cbp in (JSJ) and making use of (j3J) and (jHJ), we finally have 

<Z d-1 d d-2 

l , 

^ m (/ i i 

6=0 /3=0 



f/fd Ti) E E^ 1 + E £(p* - pm-i) (ft* - 

;. n ii h—Q 0=0 

d d-2 , 

£pm-i£ 



d d-2 

EE p ^ ( u,,i 

6=0 ,3=0 



1 



d+1" 



6=0 



)3=0 



Pb,d-\ 



d(d+ 1) 



b=0 

d d-2 

E E ; ( 

6=0 /3=0 
d d-1 

££*v (*v 

6=0 (3=0 



1 



1 



E^M-i 

6=0 



d + 1 



j-n 



6, d-1 



(16) 



This completes the proof of proposition. ■ 

Note that reconstruction formula fl9]) is not unique and many alternative expressions 
can be found, but we will use this formula in view of its symmetry Also, an immediate 
consequence of formula (Q is that if MUBs exist then any qudit state can be represented 
by a single probability distribution {(d+ l)~ 1 p aa }- From this point of view, MUB-based 
representation of qudit states is a partial case of an inverse spin-s portrait method [9] 
with an extra requirement on the symmetry. 



4. Star-product scheme 

Any operator A acting on Hilbert space of quantum states can be alternatively described 
by a symbol /a( x ) which is a function of a particular set of variables x. The relation 
between A and /a( x ) is defined through 

/4x)=Tr[i*7(x)], A= /dx/U(x)D(x), (17) 

where t/(x) and -D(x) are dequantizer and quantizer operators, respectively; an explicit 
form of the sign of integration J dx depends on the scheme used. It is worth noting that 
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substituting the second equality ( TT7|) for A in the definition of symbol of operator, we 
readily obtain that a function 2)(x!,x) = Ti[D(^xi)U(x)] has a sense of delta-function 
on symbols, i.e. / dxiD(xi, x)/ A (xi) = /a(x). 

Using such a formalism, one can deal with symbols of operators instead of operators 
themselves. The rules of addition and multiplication for operators are then transformed 
into the following rules for symbols: 

C=A+B: / c (x) = / A (x)+/ fl (x), (18) 

C = cA: / C (x) = c/ A (x), (19) 

C = AB: /c(x) = (/A*/fl)(x) s (20) 

where by star we denote a star product of symbols defined through 



(/a*/b)(x) = yy dx 1 cZx2/A(x 1 )/ B (x 2 )K(x 1 ,X2 ) x), (21) 

K(x 1; x 2 ,x) = Tr[£>(x 1 )Z)(x 2 )^(x)]. (22) 

The term i^(x 1 ,x 2 ,x) is usually referred to as star-product kernel [21] [22] . As the star 
product is associative, from the relation f A -kf B -kf c = f A * (f B * f c ) = (f A * f B ) * f c 
we readily obtain a star-product kernel of 3 symbols K^fax, x 2 , x 3 , x) as well as an 
additional requirement on the star-product kernel 

^)(x 1 ,x 2 ,X3,x)=/^(x 1 ,x 2 ,y)if(y,x 3 ,x) 

= f dyA'( Xl ,y,x)A'(x 2 ,x 3 ,y), (23) 
which is valid for all sets of variables xi, x 2 , X3, and x. 

4-1. Dual star-product scheme 

The star-product scheme of the form 

/i(x) = Tr[i£(x)], A = j dx/i(x)C/(x) (24) 

is called dual with respect to (|T7j) . Apparently, dual star-product kernel reads 
i^ d (x 1; x 2 ,x) = Tr[f/(x!)?7(x 2 )Z)(x)] and satisfies the relation fl23l) . 

4-2. Relation between tomographic schemes 

Suppose that we are given two star-product schemes: (i) U{x), D(x) and (ii) W(£), 
T?(£). A relation between the corresponding symbols is 



/i(x) = y </::(^)A:: ,($. X ). /„(£) = J dx/, (x) (x , £) , (25) 

where intertwining kernels Ka-+i(£,x) and lifi-^x, £) read 

tfii-i&x) = Tr[P(£)£>(x)], 2£-a(x,£) = Tr[L>(x)ZY(£)]. (26) 

The relation between MUB-symbols and symbols of a symmetric informationally 
complete POVM (SIC-POVM) will be considered for qubits in Section [51 
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4-3. MUB star-product scheme 

Comparing MUB-scanning procedure (GO) and reconstruction procedure (|9]) with the 
scheme (JT7|) . it is not hard to see that MUB-tomography can be treated as a star- 
product scheme with the following dequantizer and quantizer operators: 

U a a n^oaj D aa l^aa ~, '. i 

a + 1 



is a 



where x = {a, a}, a = 0, . . . , d, a = 0, . . . , d — 1 and f dx. = ^2 a=0 S a =Q, i-e. x i 
set of discrete variables and integration J dx implies summation. MUB-tomographic 
symbol of any operator A acting on rf-dimensional Hilbert space is 

f A (a,a) = Tr[AfL aa ], (28) 

d d— 1 / j \ 

^ = EEW fl ' a ) n aa - —i ) . (29) 

a=0 a=0 ^ ' 

Delta-function on MUB-tomographic symbols is 
£(a, a- b, 13) = Tr[D aa U bp ) = Tr[fl aQ % 



d+l 



1 + Sa,b (5 a ,f, ~ \ ) . (30) 



d(d+l) a '° V ' P d y 
Note that the obtained delta-function contains extra terms in addition to the Kronecker 
delta symbol 8 a j8 a> p. There is no contradiction here and it can easily be checked that 
the residual part always gives zero while summation with any MUB-symbol. 

MUB star-product kernel is expressed through MUB triple product T aayb p fil = 
Tr[n aQ n fe( gIl C7 ] as follows: 

K(a, a; b, (3; c, 7) = Tr[D aa t) bl3 U n } 

_ ji . 5 a , c + 8b,c _ 5a,c^a,7 + t^6,c^,7 __ d + 2 

~ a ° MC7 d(d+l) d+l d{d+lf [ ' 

Star-product kernel K(a, a; b, j3; c, 7) necessarily meets the condition (I23l) . from 
which we derive a new relation on MUB triple product 

d d-l 
c=0 7=0 

^(1 - 8a,b) + 5a,bb a ,f^j Q(l - + 4,«^,A 

-(1 - 5 aj l) + 5 a ,l5a : \^ ^-(1 - 5 6)fc ) + tf&.Jfe^,*^ (32) 

and find an expression which relates the four-product Tr[n aa n^n fc ^fl iA ] and triple 
product 

d d-l 

Tr[n aQ ,n^nfc^fi; A ] = T aQ , )fe/3iC7 T C7ifcx iA 

c=0 7=0 

i(l - 5 Qjb ) + 5 a>b 5 a A ( ~(1 - 5 k}l ) + 5 kJ 5 x ^ . (33) 
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It is worth mentioning that the same result can be alternatively obtained by using 
dual MUB star-product kernel of the form 

K d (a, a; b, (3; c, 7) = Ti[D aa D bf} U c 



C7J 



T, 



1 



aa,b/3,c"/ 



d+ 1 \ d 



1 



'1 - 8a,b) + S a ,b5. 



(34) 



4-4- Lie algebraic structure of MUB-POVM 

The developed MUB star-product scheme enables to reveal the Lie algebraic structure 
of MUB-projectors. In fact, following the ideas of [16], let us consider a commutator 

c = [ru,n 



sa, J-J-6/3 



Ilaall^ — n^Ilaa. Since MUB-projectors are Hermitian, we obtain 

= —C. This means that the MUB-symbol of such a commutator is purely imaginary, 
that is 



T 



T, 



fc{<-- ~ ) = Tr[Cn c7J ,.„„.,, 

where J aa ,bf3, n 
and reconstruction formula ( 1291) . we readily obtain 

d d-1 



bfi,aa,cy 
d-1 



is real and satisfies the condition E 7 =o J a a,b/3,cy 



iJaa,b/3,c"fi (35) 

0. Using this condition 



C — ^ ^ iJaa,bp,c^ 
c=0 7=0 

d d-1 



Un d+1 



n aa , n 6 



c=0 7=0 



(36) 



(37) 



The latter equation means that MUB-projectors form the Lie algebra gl(<i, C), with 
iJaa,bj3,cy being structure constants. Evidently, MUB-POVM effects {E aa } satisfy 



E aa , Eb 



aa,bf),cyE c ^. 



5. MUB star-product scheme for qubits 



MUB-projectors in 2-dimensional Hilbert space can be chosen as follows: 



n 



a=0,a=0 



^ + fr x ), 



n a= o )Q ,=i 
n a= i a= \ 



n a =i,a=o = 2 (i + ^i/)) n a= i )0= i = -(/ - fry) 
n a =2,a=o = + &z), n a=2 ,Q=i = -(/ - fr z ) 



(38) 
(39) 
(40) 



2 \- ' -^/7 Ur— i,Ll — -L 2 

where <r = (a^, a y , fr z ) is a set of Pauli operators. 

The delta-function on MUB symbols for qubits is D(a, a;b, (3) = | + 8 a> b(8 a> p — |). 
The direct calculation shows that the triple product of MUB-projectors for qubits reads 

1 + 2(8 atb 5 a) p + 8 hjC 8p n + 5 Cja 5 7jQ ) 



T, 



aa,b/3,cy 



4 



- (8 a ,b + $b,c + 8c,a) 

+ i£ a bc{8 a ,o - 8 a ,i)(8p,o - 8p,i)(8y t0 - 5 7) i) 



(41) 
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Table 1. Sign function S(k;a,a). 



k 00 01 10 11 20 21 

1 + - + - + - 

2 + - - + - + 

3 - + + - - + 

4 - + - + + - 



where e a b c is Levi-Civita symbol. Substituting the obtained triple product in formulas 
( 13T|) . (134"]) . and ( )33|) . it is easy to calculate the ordinary and dual MUB star-product 
kernels as well as the four-product for qubits. 

5.1. Relation to SIC star-product scheme 

A star-product scheme based on symmetric informationally complete POVM is 
considered in [24J. A self-dual star-product scheme that is very similar to SIC 
star-product scheme is considered in [25]. In case of qubits, SIC-projectors read 
V k = §(/ + (<x • n fc )), k = 1,...,4, where ni = ^(1,1,1), n 2 = ^=(1,-1,-1), 
n 3 = ^=(—1,1,-1), and n 4 = ^=(—1,-1,1). The dequantizer is £4 = \Ph and the 
quantizer is T>k = 37^ — I- A calculation of intertwining kernels (126]) between MUB and 
SIC star-product schemes yields 

-K'sic^mubO; a, a) = - (l + V3S(k] a, atfj , (42) 

^MUB^sic(a, a; k) = ( x + v^^A;; a, a)) , (43) 

where S(k; a, a) is a sign function taking values ±1 in accordance with table 1. We hope 
that analogues simple relations between MUB and SIC quantization schemes exist in all 
prime dimensions (MUBs and SIC-POVMs are compared also in [251 12~T]). 

6. MUB-tomography and Stern-Gerlach measurements 

On passing a beam of spin-j particles in a state p through Stern-Gerlach apparatus 
oriented along z-axis, we are able to measure probabilities to find particles in each 
splitted beam, i.e. in the state \jm), where m = is the spin projection on 

z-axis. States {\jm)y m= _j form the first basis in d- dimensional Hilbert space with 
d = 2j + 1. Suppose a magnetic field B a is applied to spin particles before they are 
passed through the Stern-Gerlach magnetic field gradient. This results in a unitary 
transformation of the initial state p — > u\pu a . The probabilities of outcomes read 
p a {fn) = (jm\ulpu a \jm) . On the other hand, p a (m) = p aa = (aa\p\aa), where \aa) = 
&o|j, m — a ~ j)i a — 0, • • • , 2j. States {\aa}}^ =0 form a new basis in Hilbert space, 
with parameter a being a label of this basis. Thus, applying different magnetic fields B a , 
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a — 0, . . . , 2j + 1, we obtain a set of 2j + 2 bases {\act)}J =0 . If unitary transformations 
u a satisfy an additional condition \ {jm\u\ui,\jm')\ = (2j + 1) _1 (1 — 5 a ,b) + S a ,bSm,m', 
then the constructed bases are MUBs. 

However, application of constant magnetic fields {Bj^j, 1 gives rise to unitary 
transformations {ia}^ 1 of the group SU(2). Conversely, the MUB-condition (CQ) can 
be only met when u a G SU(N) with N = 2j + 1. As a result, MUBs can be constructed 
via the conventional Stern-Gerlach technique for qubits (j = |) only For higher spins, 
the minimal necessary number of unitary transformations u a G SU{2) is known to be 
4j + 1 [9j which is greater than 2j + 2. 

7. Conclusions 

To conclude we summarize the main results of the paper. 

Starting from peculiarities of mutually unbiased bases, we show that whenever 
MUBs exist they can be used in quantum state tomography Then we develop 
the MUB-tomographic-probability representation of quantum states by considering 
MUB-projectors within the framework of a star-product scheme. Dequantizer and 
quantizer of MUB star-product scheme are shown to be easily expressed through MUB- 
projectors. This takes place due to a high symmetry of MUBs. For the constructed 
MUB quantization scheme, ordinary and dual star-product kernels are calculated and 
expressed through triple product of MUB-projectors. Employing the specific rules of 
the star-product kernel, we derive a new relation on triple- and four-products. Applying 
the MUB star product scheme to a commutator of MUB-projectors, we reveal the 
Lie algebraic structure of MUB-projectors and find structure constants. The obtained 
results can be used both in seeking and classification of MUBs in higher dimensions. 

Example of qubits is considered in detail. In particular, an explicit form of all 
star-product characteristics is obtained. Intertwining kernels between MUB and SIC 
star-product schemes are found. These kernels can be used in order to find SIC- 
POVMs whenever MUBs exist, for instance, in all prime dimensions. Finally, an 
implementation of MUB-tomography via Stern-Gerlach apparatus is discussed. The 
conventional experiment turns out to be appropriate for MUB-tomography of qubits. 
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